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I. 


Introduction 


Wind  measurements  in  the  mesosphere  have  been  far  less  numerous  than 
in  the  stratosphere  below  or  in  the  lower  thermosphere  above.  Rawinsondes 
and  small  meteorological  rockets  provide  relatively  detailed  wind  data  below 
about  65  km,  while  meteor-wind  radars  routinely  provide  data  from  about 
90-105  km  at  several  stations  globally.  As  a result,  knowledge  of  the 
circulation  patterns  from  65-90  km  has  been  largely  based  on  a relative 
handful  of  rocket  grenade  experiments  (CIRA,  1972).  The  wind  measurements 
made  at  Saskatoon  (52°N,  106°W)  by  the  radio  reflection  technique  (Gregory 
and  Stephenson,  1972),  extending  from  about  65-118  km,  are  thus  of  special 
interest  as  they  help  connect  knowledge  of  the  wind  fields  of  the 
stratosphere  and  lower  thermosphere. 

The  ground-based  radio  reflection  technique  was  developed  many  years 
ago,  but  only  recently  has  it  been  practical  to  collect  and  process  large 
quantities  of  wind  data.  The  automated  data  collection  system  at 
Saskatoon  has  been  refined  so  that  soundings  with  3 km  height  resolution 
from  65-118  km  can  now  be  obtained  at  5-minute  intervals  on  a nearly 
continuous  basis.  Processing  of  the  raw  data  to  calibrate  them  as  winds 
remains  a formidable  task,  and  with  the  increased  number  of  soundings  a 
backlog  of  unprocessed  data  accumulated. 

This  report  describes  an  effort  to  process  the  backlog  of  raw  data 
and  to  analyze  all  available  wind  data  from  Saskatoon.  The  data  processing 
effort  is  documented  in  detail  in  Section  II,  while  the  analysis  results  are 
given  in  Section  III.  The  physical  principles  of  the  measurements  and 
restrictions  to  the  data  are  discussed  by  Gregory  and  Manson  (1975).  It  is 
noted  here  only  that  although  the  quoted  altitudes  above  about  105  km  are 
too  large  by  up  to  5 km  due  to  the  existence  of  group  delay  effects,  data 
are  used  here  as  reported. 
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A copy  of  Che  wind  calculation  program  was  supplied  by  D.  G. 
Stephenson  of  the  University  of  Saskatchewan,  Saskatoon.  The  program  had 
been  translated  from  IBM  FORTRAN  and  run  at  the  National  Center  for 
Atmospheric  Research  (NCAR)  on  a CDC  7600.  The  program  was  examined  to 
learn  the  wind  finding  procedure,  and  timing  runs  were  made  to  determine 
any  parts  of  the  program  that  could  be  rewritten  to  significantly  decrease 
the  execution  time.  The  timing  runs  were  made  with  the  data  from  10  June 
1975  at  1800Z.  This  seemed  to  be  a typical  sounding  with  winds  found  at 
seven  of  the  possible  23  levels. 


The  first  timings  were  made  using  the  CDC  RUN  compiler.  This  compiler 
was  needed  to  interface  with  the  input  and  unpacking  routines  supplied  by 
NCAR.  Two  parts  of  the  program  were  identified  as  potential  candidates 
for  revision:  correlation  and  unpacking.  Table  1 shews  that  of  a total 

execution  time  of  34  seconds,  27  seconds  were  spent  in  computing  correla- 
tions and  3.7  seconds  were  spent  unpacking  the  input  data.  Execution  times 
given  here  are  CDC  6600  central  processor  seconds. 

The  correlation  section  was  the  most  obvious  candidate  for  revision. 
This  section  consisted  of  four  nested  loops  of  which  the  innermost  computed 
the  sums  of  lagged  products.  The  innermost  loop  was  rewritten  to  fit  into 
the  6600  instruction  stack  which  is  a group  of  registers  that  hold  instruction 
words.  Loops  that  fit  instack  can  execute  quickly  since  the  instruction 
words  need  to  be  read  from  memory  only  once.  This  reformulation  cut  the 
correlation  time  from  27  to  12  seconds. 


At  this  point,  the  program  was  modified  slightly  so  that  it  could  be 
compiled  with  FORTRAN  EXTENDED,  a compiler  that  is  capable  of  extensive 
code  optimization.  Now  the  total  execution  time  was  down  to  10  seconds  of 
which  3 were  for  correlations  and  2.6  were  for  unpacking. 
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The  unpacking  portion  of  the  program  was  modified  by  replacing  the 
calls  to  the  NCAR  supplied  general  unpacking  routine  with  instructions 
tailored  to  this  particular  unpacking  problem.  With  these  instructions, 
which  use  the  hardware  shifts  and  boolean  operations,  the  unpacking  time 
was  cut  from  2.6  to  0.6  seconds. 

After  the  program  optimization  had  been  completed,  a program  incor- 
porating a different  wind  finding  method  was  supplied  by  Dr.  Stephenson. 

An  examination  of  the  new  program  revealed  that  only  the  unpacking  section 
could  readily  be  made  more  efficient.  The  changes  that  had  been  made  to  the 
unpacking  logic  of  the  old  program  were  applied  to  the  new  program  also. 

The  total  execution  time  of  the  new  program  was  now  4.8  seconds. 

The  26  tapes  listed  in  Table  2 were  processed  to  retrieve  winds. 

Each  file  contained  the  raw  data  for  one  sounding.  The  first  step  in  the 
processing  was  to  reformat  the  nine-track  tapes  sent  from  Saskatoon. 

Several  of  the  shorter  input  tapes  were  consolidated  onto  one  reel  for  the 
winds  processing.  The  d were  not  quality  checked  at  this  point  but 
were  merely  copied  t nat  suitable  for  the  CDC  6600  on  which  the  winds 

would  be  computed 

The  reformatted  tapes  were  used  as  input  to  the  winds  program  which 
had  been  modified  to  make  it  more  efficient.  The  outputs  of  the  program 
were  cards  and  listings  which  were  actually  saved  on  tape  rather  than  being 
punched  or  printed.  One  card  contained  the  wind  for  one  height  level  along 
with  identification  and  various  intermediate  parameters  and  flags.  A card 
was  generated  sometimes  even  if  the  program  could  not  carry  the  computations 
far  enough  to  find  a wind  value.  The  listings  contained  the  winds  along 
with  the  results  of  many  intermediate  calculations.  If  the  program  was 
unable  to  produce  a wind  for  a level,  the  listing  contained  an  explanation 
of  the  problem  with  the  data. 

The  final  step  in  the  winds  processing  was  to  consolidate  the  output. 
The  cards  were  sorted  chronologically  and  rewritten  to  a blocked  nine  track 
tape  to  be  sent  to  Saskatoon.  A total  of  33,192  cards  were  produced  from 
the  tapes  listed  in  Table  2.  Dr.  Stephenson  also  sent  the  October  to 
December  1975  cards  to  Oxford  University  to  be  compared  with  data  from  the 
Pressure  Modulator  Radiometer.  The  tapes  containing  the  listabie  output 


> 4**.:  ...  «. 


• * 


At 


* t .V  . . . ■ 


4 


were  also  rewritten  in  a format  suitable  for  use  at  the  University  of 
Saskatchewan.  These  tapes  contained  a total  of  1,905,076  printed  lines  of 
output . 

III.  Analysis 


Based  on  the  frequency  of  observations  and  the  time  periods  over 
which  they  were  taken,  the  analysis  of  Saskatoon  wind  variability  conveniently 
is  divided  into  three  classes:  variations  longer  than  a month,  variations 

longer  than  a day,  and  gravity/ tidal  wave  variations.  Before  discussing 
the  results,  the  available  data  summarized  in  Table  3 are  discussed. 

Observations  during  1969-1973  were  taken  from  publications  (e.g., 
Manson,  et  al.,  1974),  research  reports  (e.g.,  Gregory,  et  al. , 1973),  or 
from  working  charts  furnished  by  A.  Manson.  Whenever  there  was  data  for 
more  than  one  hour  each  day,  only  noon  observations  were  used.  The  plotted 
numbers  from  the  reports  were  transcribed  directly  to  punch  cards  along 
with  an  estimated  date  for  each  given  vertical  profile.  In  Table  3,  the 
number  of  hourly  medians  indicated  for  1969-1973  is  the  number  of  values 
extracted  from  the  sources  mentioned.  Because  the  number  of  soundings  on 
which  each  value  was  based  was  not  always  available,  each  value  was  counted 
as  one  observation  whether  presented  as  a daily  or  weekly  median. 

All  other  data  reflected  in  Table  3 were  available  as  individual 
soundings.  For  the  period  January  1974  through  January  1976,  daily  median 
values  at  local  noon  (18-19  GMT)  only  were  used.  The  hourly  medians  were 
formed  by  taking  components  of  the  individual  wind  vectors  at  each  level, 
and  then  selecting  the  median  value  of  this  sample  for  each  component. 

Although  the  use  of  hourly  medians,  as  opposed  to  means,  reduces  the 
occurrence  of  wild  values,  further  quality  control  seemed  necessary.  The 
deviations  from  the  annual  and  other  long  period  variations  were  obtained 
by  high-pass  filtering  the  data  at  each  level.  The  high-pass  filter  had 
507o  response  at  period  61  days.  If  a filtered  value  was  more  than  three 
standard  deviations  removed  from  the  mean  of  the  filtered  values,  then 
the  corresponding  data  point  in  the  original  series  was  discarded.  As  the 
standard  deviation  of  the  filtered  values  varied  by  less  than  10%  between 
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summer  and  winter  at  most  levels,  the  quality  control  procedure  was  applied 
to  the  entire  period,  January  1974  to  January  1976,  without  regard  for 
season.  About  4%  of  the  observations  at  each  level  were  thus  discarded. 

Profiles  were  taken  at  5-minute  intervals  throughout  the  10-day 
intensive  runs  in  August  and  October,  1976,  although  at  least  50%  of  the 
total  possible  observations  are  missing  from  each  series  (Table  3). 

Seasonal  or  long-term  synoptic  variations  were  removed  from  each  intensive 
observation  run  by  subtracting  a parabolic  trend  line  fit  to  the  data  by 
the  least-squares  method.  These  data  were  quality  controlled  by  the  same 
procedure  as  above,  excap t that  the  filter  had  50%  response  at  period  4.5 
hours . 


A.  Long-period  oscillations 

This  section  covers  the  first  four  harmonics  of  the  annual  wave  and 
the  quasibiennial  oscillation  (QBO) . Some  of  these  oscillations  can  be 
seen  in  the  time  series  plots  of  the  observed  wind  speed  presented  next, 
and  the  amplitudes  and  phases  of  the  waves  from  periodic  analysis  are  given 
later. 


Monthly  mean  profiles  of  the  zonal  wind  speed  (u)  are  shown  in 
Figure  1.  There  are  no  noteworthy  changes  from  the  results  of  Gregory,  et  al. 
(1975),  which  are  based  on  data  from  1969-1973,  although  the  two  additional 
years  of  data  used  here  doubles  the  number  of  hourly  medians  at  most  levels. 

The  annual  wave  is  clearly  seen  in  Figure  1 and  has  a phase  reversal  near 

100  km  with  winter  westerlies  below  and  summer  westerlies  above  that  level. 

In  the  overlap  region  of  these  two  annual  waves,  i.e.,  from  90-100  km,  a 
semiannual  wave  is  present  as  there  are  four  reversals  per  year.  The  semi- 
annual and  shorter  period  waves  are  not  easily  seen  by  inspection  at  other 

levels . 


It  is  interesting  that  a jet  of  winter  westerlies  is  not  found  in 
the  lower  mesosphere  in  Figure  1.  As  the  mean  speed  in  the  jet  is  over  60 
ms  ^ both  from  rocketsonde  data  (Belmont,  et  al . , 1975)  and  grenade  data 
(Theon,  et  al.  , 1973),  the  relatively  small  winter  speeds  below  80  km  in 
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Figure  1 may  reflect  an  instrument  bias  toward  low  speeds.  Another  possible 
reason  that  a winter  jet  is  lacking  is  that  only  noon  data  are  used  in 
Figure  1,  while  the  rocket  and  grenade  data  at  Churchill  (59°N)  are  taken 
at  all  hours  of  the  day  (Nastrom  and  Belmont,  1975;  Theon,  et  al.,  1973). 

Thus,  the  diurnal  tide  is  included  in  the  noon  data  used  here  and  seasonal 
variations  of  amplitude  or  phase  of  this  tide  may  bias  the  monthly  mean 
noon  wind  speeds  such  that  the  annual  wave  below  about  80  km  appears  too 
small.  There  is  not  sufficient  data  at  this  time  to  verify  possible 
seasonal  changes  in  the  diurnal  tide  below  80  km,  and  the  large  changes 
above  80  km  described  by  Fellous,  et  al.  (1975)  should  not  be  extrapolated 
downward.  Below  80  km  the  diurnal  tide  has  relatively  large  amplitude 
seen  in  the  height-time  section  of  mean  hourly  data  for  October  13-22,  1976, 
in  Figure  2.  Seasonal  changes  in  phase  or  amplitude  would  clearly  affect 
the  annual  cycle  of  noon  winds  significantly.  As  a thorough  analysis  of 
tidal  winds  during  each  season  is  being  made  by  the  Saskatoon  group,  tides 
will  not  be  discussed  further  here.  Until  their  results  are  available,  it 
is  possible  only  to  acknowledge  that  tidal  biasing  may  be  important. 

Finally,  possible  biasing  by  the  diurnal  and  semidiurnal  tides  must 
be  taken  into  account  if  the  present  data  are  used  for  developing  models 
of  the  general  circulation.  Otherwise  the  resulting  models  may  be  dynami- 
cally unstable  (Schoeberl  and  Zalesak,  1976). 

The  mean  monthly  noon  meridional  winds  (v)  are  given  in  Figure  3. 

The  basic  pattern  is  similar  to  that  of  Gregory,  et  al.  (1975).  These  results 
are  included  for  completeness  and  will  not  be  discussed  further. 

A QBO  in  zonal  wind  speed  in  the  lower  thermosphere  has  been  noted 
by  Sprenger,  et  al.  (1975)  and  in  the  airglow  of  the  lower  thermosphere  by 
Fukuyama  (1977).  Long-period  oscillations  are  also  found  at  Saskatoon 
(Figure  4).  Twelve-month  running  means  of  monthly  mean  noon  winds  are 
plotted  in  Figure  4 and  the  periods  in  months  between  successive  maxima  or 
minima  are  noted.  Missing  data  cause  the  filter  to  be  quite  leaky,  and 
the  results  in  Figure  4 can  thus  only  be  used  qualitatively.  An  interpolation 


technique  perhaps  could  have  been  used  to  obtain  a complete  series,  and  a 
better  filter  could  then  have  been  used  to  more  precisely  isolate  the  QBO. 
However,  uncertainties  then  would  merely  have  been  imbedded  in  the  interpo- 
lation process  which  could  be  no  more  reliable  than  the  original  data,  so 
the  present  method  was  deemed  adequate  for  the  present  purpose.  The  point 
here  is  that  variations  with  a QBO  period  are  present  at  Saskatoon.  The 
QBO  in  thermospheric  wind  is  not  confined  to  northern  mid-latitudes,  for  a 
QBO  is  also  apparent  at  100  km  in  the  12-month  running  means  of  meteor 
wind  data  at  Adelaide  (35°S)  (Figure  5)  from  Elford  (1977).  The  Adelaide 
winds  are  based  on  approximately  once-weekly  soundings. 

In  order  to  estimate  the  amplitudes  and  phases  (time  of  maximum 
eastward  speed)  of  the  QBO  and  of  the  first  four  harmonics  of  the  annual 
wave,  periodic  analysis  of  monthly  mean  zonal  wind  speed  was  made  using 
the  technique  of  Belmont  and  Dartt  (1973).  Tests  with  sinusoidal  waves 
having  periods  ranging  from  20  to  35  months  showed  that  a 26  month  period 
most  reduced  the  variance  and  gave  the  smallest  RMS  amplitude  errors  for 
this  period  of  record.  Vertical  profiles  of  the  amplitude  and  phase  of  the 
26-month  QBO  are  given  in  Figure  6a.  Below  118  km,  QBO  amplitudes  range 
from  1 to  7 ms  \ in  good  agreement  with  Sprenger,  et  al.  (1975)  who  report 
amplitudes  of  about  3 ms  at  90-95  km  near  53°N.  The  phase  of  the  QBO  is 
highly  variable,  and  perhaps  the  only  reliable  feature  is  a 180°  shift  between 
94  and  97  km. 

The  amplitude  of  the  annual  wave  (Figure  6b)  has  a pronounced  minimum 
near  100  km.  Below  100  km  the  phase  date  is  near  the  winter  solstice, 
and  above  100  km  it  is  near  the  summer  solstice.  Above  105  km  the  present 
amplitudes  are  slightly  larger,  and  the  phase  dates  about  a month  later  than 
the  corresponding  results  of  Groves  (1972).  The  unexpectedly  small 
amplitudes  below  about  80  km  were  discussed  in  connection  with  Figure  1. 
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The  amplitude  of  the  semiannual  wave  (Figure  6c)  has  three  maxima, 
in  agreement  with  Groves  (1972).  The  phase  dates  are  also  similar  to  those 
of  Groves,  but  the  very  regular  pattern  of  propagation  with  height  is  much 
more  evident  in  the  present  display  and  does  not  seem  to  have  been  noted 
previously.  The  straight  line  on  the  phase  diagram  (fitted  by  eye)  has  a 
slope  of  8 km  month  If  downward  phase  propagation  implies  upward  energy 

propagation,  then  the  region  of  excitation  of  this  semiannual  wave  is  near 
70  km.  Belmont,  et  al.  (1973),  have  hypothesized  that  particle  precipitation 
effects  in  the  auroral  zone  influence  the  ozone,  and  hence,  thermal  and  wind 
fields  on  a semiannual  basis.  Recent  observational  (Randhawa,  1976;  Heath, 
et  al.,  1976)  and  theoretical  (Thorne,  1977)  papers  lend  support  to  the 
particle  precipitation  hypothesis,  and  the  present  results  could  be  explained 
by  that  hypothesis.  However,  there  could  be  other  explanations  for  this 
semiannual  wave.  For  example,  the  wind  reversals  associated  with  mid-winter 
sudden  warmings  could  induce  a semiannual  wave  with  equinoctial  phase  dates, 
as  found  in  Figure  6c  near  70  km.  However,  sudden  warmings  should  induce 
nearly  the  same  phase  at  all  levels. 

The  3-month  wave  can  be  identified  as  the  second  harmonic  of  the 
semiannual  wave.  Although  its  amplitude  (Figure  6d)  is  less  than  5 ms  ^ 
at  most  levels,  its  phase  progression  with  height  appears  exactly  double 
that  of  the  semiannual  wave.  This  is  apparently  the  first  time  an 
organized  3-month  wave  has  been  detected  in  the  mesosphere  or  lower 
the rmosphere. 

The  terannual  wave  has  amplitudes  less  than  5 ms  ^ at  most  levels 
and  has  no  clear  pattern  of  phase  progression  with  height,  so  is  not  shown 
here. 


B.  Planetary  waves 

Planetary  waves  are  those  oscillations  with  periods  from  2 days  to 
over  a month.  At  middle  latitudes  in  the  lower  atmosphere  these  waves  are 


generally  identified  as  Rossby  waves.  At  higher  levels,  peaks  in  the 
variance  power  spectrum  of  meteor  wind  data  have  been  found  at  about  5 
and  10  days  (Teptin,  1972).  Stimulated  by  the  modeling  results  of  Geisler 
and  Dickinson  (1976),  Fraser  (1977)  has  examined  the  5-day  wave  in  iono- 
spheric absorption  in  detail.  Muller  and  Kingsley  (1974)  and  Glass,  et  al. 
(1975)  have  also  studied  the  5-day  wave  in  meteor  winds.  A 5-day  wave  could 
thus  be  expected  in  the  present  wind  data. 

Another  oscillation  which  might  be  anticipated  in  the  lower 
thermosphere  is  the  approximately  27-day  period  of  the  solar  rotation. 

This  is  suggested  by  the  well-known  27-day  cycle  in  geomagnetic  activity  and 
the  report  by  Banks  (1977)  that  significant  heating  may  occur  in  the  lower 
thermosphere  of  the  auroral  zone  during  disturbed  periods. 

To  examine  whether  these  or  any  other  possible  waves  are  in  the 
Saskatoon  wind  data  variance  power  spectra  are  given  in  Figures  7 and  8. 

These  results  are  based  on  local  noon  daily  data  averaged  over  12  km  height 
intervals,  although  spectra  were  also  taken  at  all  individual  levels,  as 
discussed  later.  Because  of  possible  seasonal  differences,  the  analysis  was 
made  for  both  summer  (May-September)  and  winter  (October-April) . The 
statistical  significance  of  spectral  peaks  can  be  tested  by  comparison  with 
the  95%  confidence  limit  of  a first-order  Markov  "null"  model  shown  by  dashed 
lines  in  Figures  7 and  8.  The  assessment  of  significar.ce  is  thus  independent 
of  the  coordinates  used  for  presentation.  To  help  account  for  missing 
observations,  the  number  of  data  points  was  taken  to  be  the  number  of  pairs 
at  lag  one  day.  The  maximum  lag  used  in  forming  the  autocorrelation  function 
was  61  days. 

During  winter  (Figure  7),  marginally  significant  peaks  at  period 
4.8  days  are  found  in  both  u and  v at  85-97  km.  A nearly  significant  peak 
is  found  in  v at  the  same  period  at  70-82  km  and,  from  cross-spectral 
analysis,  has  a phase  lag  of  -3°  relative  to  85-97  km  although  the  coherency- 
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square  is  only  0.15.  No  peaks  are  found  in  the  spectral  band  centered  near 
the  solar  rotation  period  of  27  days.  Between  2 and  3 days  period  several 
peaks  in  Figure  7 exceed  the  95%  level,  but  the  variety  of  periods  at  which 
they  occur  does  not  permit  associating  them  with  the  2-day  wave  discussed 
by  Glass,  et  al.  (1975).  At  70-82  km,  peaks  are  found  in  both  v and  u at 
period  15.25  days,  although  only  the  former  is  significant. 

The  model  results  of  Geisler  and  Dickinson  (1976)  suggest  that  the 
largest  5-day  wave  at  this  latitude  (52°N)  should  be  found  in  the  v-component 
in  summer  above  80  km.  In  Figure  8,  a nearly  significant  peak  is  found  at 
period  4.9  days  at  85-97  km,  and  peaks  at  4.4  and  7.2  days  are  seen  at 
100-112  km.  In  the  spectral  band  centered  at  30.5  days,  peaks  are  found  at 
all  levels  in  v and  at  85-97  km  in  u,  although  none  exceeds  the  95%  level. 
Significant  peaks  at  periods  slightly  over  2 days  are  found  at  all  levels  in 
both  components,  and  probably  correspond  to  the  2 -day  wave  of  Glass,  et  al 
(1975). 

Waves  which  have  rapid  phase  changes  with  height  or  which  exist  over 
only  a portion  of  one  of  the  height  intervals  used  above  might  be  obscured 
due  to  averaging.  To  examine  that  possibility  the  spectral  analysis  results 
for  summer  at  individual  levels  are  summarized  in  Table  4.  A dot  is  entered 
in  the  table  if  a spectral  peak  is  found  at  a given  level  and  period,  while 
an  S is  entered  if  the  peak  exceeds  the  95%  level.  Forcing  at  a given  period 
may  be  expected  to  induce  a spectral  peak  at  several  adjacent  height  levels, 
although  all  peaks  may  not  necessarily  be  significant.  Possible  frequency 
shifts  with  height  are  ignored  here.  Using  the  ad  hoc  criterion  that  a peak 
at  fi'-?  adjacent  levels  may  suggest  some  organized  forcing,  then  only  the  30.5 
day  period  and  periods  near  4.5,  3,  and  2 days  would  be  judged  physically 
meaningful.  Corresponding  results  for  winter  only  confirm  the  findings 
discussed  before,  so  are  not  given  here. 


. ! "Z?  , »i'  ; — - 


. «v 


11 


C.  Short  period  motions 


* J 


This  section  considers  those  motions  with  periods  less  than  a few 
hours,  thus  including  gravity  waves  and  turbulence.  Tides  were  discussed 
briefly  in  Section  IIIA  and  are  not  considered  further  here,  although  it 
must  be  recognized  that  higher  harmonics  of  the  diurnal  and  semidiurnal 
tides  could  be  included  in  the  present  results.  One  aim  of  this  study 
is  to  determine  if  a sampling  interval  longer  than  five  minutes  could  be 
used  to  adequately  resolve  the  motion  systems  of  the  mesosphere  and  lower 
thermosphere.  If  there  is  a distinct  gap  in  the  spectrum,  then  observations 
could  be  spaced  at  twice  the  frequency  of  the  gap. 

In  order  to  examine  the  partitioning  of  the  variance  among  motions 
with  periods  from  10  minutes  to  6 hours,  spectra  of  the  entire  10-day 
series  of  individual  5-minute  observations  during  August  and  October,  1976, 
were  computed.  Inspection  of  those  results  indicated  that  prewhitening 
should  be  applied.  Straightforward  removal  of  the  mean  daily  variation 
(e.g.,  see  Figure  2)  was  not  used  because  of  the  well-known  day-to-day 
variation  of  the  tides.  Also,  a high-pass  filter  was  not  used  because 
differing  amounts  of  missing  data  under  the  filter  induce  differing  filter 
response  characteristics  which  could  hamper  the  interpretation  of  any 
results.  Thus,  at  each  level  each  10-day  series  was  divided  into  6-hour 
segments.  Segments  with  less  than  507o  complete  data  were  discarded,  and 
each  of  the  remaining  segments  was  detrended  by  regression  with  a parabolic 
trend  line  in  an  effort  to  remove  tidal  and  synoptic  variations.  An  average 
of  19  segments  were  used  at  each  level.  Lagged  covariances  up  to  maximum 
lag  three  hours  ware  determined  for  each  segment,  and  then  averaged  over  all 
segments.  The  variance  power  spectra  of  these  stacked  covariances  are 
presented  in  Figure  9 and  10  for  the  94  and  106  km  levels,  which  are  typical 
examp les. 


Vertical  lines  are  entered  in  Figures  9 and  10  at  the  periods  where 
507o  of  the  total  variance  is  on  either  side  of  the  line.  The  lines  are 
always  found  between  24  and  34  minutes  period.  Thus,  although  gravity 
waves  with  periods  one  hour  or  greater  could  be  resolved  with  observations 
spaced,  say,  10  or  20  minutes  apart,  any  estimates  would  be  seriously 
aliased  by  the  relatively  large  variance  contained  in  the  period  range 
10-25  minutes.  Possible  aliasing  of  the  present  results  by  motions  with 
less  than  10  minutes  period  cannot  be  assessed  directly,  but  must  be 
expected. 

The  variations  with  periods  greater  than  about  half  an  hour  have  been 
identifiedas  gravity  waves  by  Manson,  et  al.  (1974),  although  Teptin  (1975) 
prefers  to  describe  these  variations  as  turbulence.  Gravity  waves,  although 
individually  organized,  are  expected  to  be  sporadic  and  polychromatic. 

Thus,  if  a long  period  of  data  is  analyzed,  as  used  initially  here  or  as 
Teptin  used,  any  gravity  wave  influence  would  be  smeared  out,  leaving  the 
impression  that  only  turbulent  motions  exist.  For  example,  the  average 
correlation  coefficient  at  lag  one  (i.e.  , five  minutes)  at  94  km  during 
the  October,  1976  period  is  0.28  for  u and  0.08  for  v.  If  the  only  variance 
were  that  due  to  a sinusoid  of  period  one  hour , the  lag  one  correlation 
coefficient  would  be  0.866;  and  if  the  period  were  one-half  hour  the  lag 
one  correlation  coefficient  would  be  0.50.  The  observed  average  coefficients 
are  too  small  to  attribute  solely  to  a single  wave  action.  However,  the  lag 
one  correlation  coefficients  over  individual  six-hour  segments  (Table  5) 
are  occasionally  large  enough  to  be  due  to  waves.  On  other  occasions  they 
are  negative,  indicating  very  high  frequency  (turbulent)  motions.  The 
point  is  that  while  statistics  of  large  samples  are  usually  desirable, 
multiple  and  sporadic  phenomena  can  be  superimposed  causing  ambiguity. 
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IV.  Summary  and  Conclusions 

1.  A weak  quasibiennial  oscillation  of  less  than  7 ms  ^ from  65-115  km 
agrees  with  independent  estimates  at  other  locations.  The  phase  shows  a 
180°  shift  near  95  km. 

2.  The  annual  wave  is  about  20  ms  ^ up  to  88  km  and  above  105  km  with 
a pronounced  minimum  near  100  km.  Below  100  km  the  maximum  occurs  at  the 
winter  solstice  and  above  100  km  near  the  summer  solstice.  The  data  and 
the  results  given  here  may  be  biased  by  the  seasonal  variation  in  the 
diurnal  tidal  component,  as  data  are  taken  near  local  noon.  Tidal  analyses 
are  being  made  elsewhere. 

3.  The  semiannual  wave  has  three  maxima  near  10  ms  ^ at  70,  97,  and  118 
km,  agreeing  with  Grove's  model.  The  phase  displays  a remarkably  consistent 
downward  progression  of  8 km  month  ^ from  118  km  in  October  to  70  km  in 
April.  These  dates  suggest  either  influence  of  auroral  particle  precipitation 
maxima,  and  consequent  effects  on  ozone,  or  the  influence  of  heating  at 

high  latitudes  by  sudden  warmings  in  winter,  or  both. 

4.  A weak  3-month  wave  of  less  than  7 ms  ^ has  double  the 
phase  descent  rate  of  the  semiannual  wave,  of  which  it  is  the  second 
harmonic.  This  appears  to  be  the  first  notice  of  a 3-month  wave. 

5.  No  consistent  terannual  wave  phase  was  observed  and  amplitudes  were 
less  than  5 ms  \ 

6.  Among  planetary  scale  waves,  marginally  significant  peaks  at  4.8 
days  were  found  in  both  u and  v at  85-97  km  in  winter,  and  at  4.4  days  in 

v at  70-82  km.  In  summer  an  approximately  five-day  wave  appears  at  85-97  km 
and  100-112  km.  In  addition,  in  simmer  a spectral  peak  at  7.2  days  at 
100-112  km,  is  found.  Spectral  peaks  near  27  days  period  are  not  found  in 
winter,  but  small  spectral  peaks  are  found  in  summer  at  all  levels  for  v 
and  at  middle  levels  for  u.  Peaks  near  two  days  period  are  seen  at  all 
levels  in  u and  v in  both  seasons. 


<1  IV  . «■*  - 


7.  From  5-minute  data,  it  is  concluded  that  possible  aliasing  by  the 
large  variance  in  the  period  range  of  10-25  minutes  requires  observations 
at  highest  possible  frequencies. 
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Table  1 


Program  execution  times  given 

in  CDC  6600  central  processor 

seconds 

Total 

Correlations 

Unpacking 

Execution 

Initial  timing 

27 

3.7 

34 

Rewrite  correlation  section 

12 

3.7 

19 

Change  to  FORTRAN  EXTENDED 

3 

2.6 

10 

Rewrite  unpacking  section 

3 

0.6 

8 

New  Program  (correlations 

- 

0.6 

4.8 

are  done  differently) 

Table  2 

Raw  data  tapes  received  from  Saskatoon  and  processed  by  CDC  (See  Table  3 
for  other  available  data). 


'•  ipe  Dates  No.  of  Files 


TIDE51 

4 

Apr 

- 

26 

May 

1974 

266 

TIDE52 

30 

May 

- 

30 

Jun 

1974 

189 

TIDE53 

4 

Jul 

- 

30 

Jul 

1974 

233 

TIDE54 

1 

Aug 

- 

23 

Aug 

1974 

245 

PARK61 

8 

May 

- 

9 

May 

1975 

24 

PARK62 

10 

May 

- 

12 

May 

1975 

36 

PARK63 

13 

May 

- 

15 

May 

1975 

37 

PARK65 

20 

May 

- 

23 

May 

1975 

24 

PARK66 

24 

May 

- 

29 

May 

1975 

83 

PARK69 

7 

Jun 

- 

11 

Jun 

1975 

60 

PARK70 

13 

Jun 

- 

17 

Jun 

1975 

61 

PARK71 

18 

Jun 

- 

24 

Jun 

1975 

72 

PARK77 

10 

Jul 

- 

16 

Jul 

1975 

88 

PARK78 

17 

Jul 

- 

22 

Jul 

1975 

107 

PARK79 

23 

Jul 

- 

29 

Jul 

1975 

72 

PARK80 

30 

Jul 

- 

8 

Aug 

1975 

121 

PARK84 

27 

Aug 

- 

5 

Sep 

1975 

141 

PARKS  5 

6 

Sep 

- 

16 

Sep 

1975 

135 

PARK86 

17 

Sep 

- 

29 

Sep 

1975 

99 

PARK87 

30 

Sep 

- 

10 

Oct 

1975 

156 

PARK88 

11 

Oct 

- 

16 

Oct 

1975 

72 

PARK89 

17 

Oct 

- 

31 

Oct 

1975 

188 

PARK90 

1 

Nov 

- 

14 

Nov 

1975 

169 

PARK9 1 

15 

Nov 

- 

24 

Nov 

1975 

144 

PARK94 

11 

Dec 

- 

30 

Dec 

1975 

203 

WND008 

13 

Jan 

- 

25 

Jan 

1976 

201 

Table  3 


Summary  of  available  data 


73  88  103  118  km 

Dates  Observational  frequency Number  of  hourly  medians 


Feb  1969  - Aug  1973 

Hourly  median  at  local  noon 

68 

102 

43 

0 

once  or  twice  each  week 

Jan  1974  - Jan  1976 

Hourly  median  at  local  noon 

324 

403 

350 

329 

each  day 

Number 

of 

observations 

13  - 22  Oct  1976  / 

Continuous  data  at  5-minute 

382 

983 

848 

491 

5 - 14  Aug  1976  | 

intervals  24  hours  per  day 

223 

145 
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Table  5 


r 


Lag  one  correlation  coefficients  at  94  km 
October  1976 


Segment 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 


u 

V 

.26 

.09 

.34 

-.03 

.19 

.04 

.28 

.64 

.35 

.32 

.04 

-.10 

.43 

-.18 

.36 

-.47 

-.07 

-.66 

.27 

.29 

.44 

.40 

-.11 

.35 

.54 

.12 

0.28 

0.08 
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Figure  2.  Consolidated  hourly  mean  zonal  wind  speeds  (ms  ) during  13-22 
October,  1976.  Winds  from  the  west  are  hatched. 
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Twelve-month  running  means  of  zonal  wind  speed  at  88  and  100 
km.  Every  other  value  is  plotted,  although  some  data  are 
missing.  As  missing  data  degrade  the  filter  response  (e.g.  , 
there  is  an  apparent  annual  wave  leakage  at  100  km  in  1971), 
these  results  must  be  used  only  qualitatively.  The  number  of 
months  between  successive  like-extrema  have  been  entered. 
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Figure  5. 


As  in  Figure  4,  except  for  Adelaide  (35°S)  at  84  and  106  km, 
The  time  series  at  Adelaide  are  complete. 
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Figure  6.  Periodic  analysis  results  at  Saskatoon,  based  on  monthly  mean 

local  noon  data  1969-January  1976;  (a)  26-month  QBO,  (b)  annual 
wave,  (c)  semiannual  wave,  (d)  3-month  wave.  The  amplitude 
(ms-*)  and  phase  (degrees)  at  each  level  are  given  along  with 
the  RMS  errors  of  estimate.  See  text. 
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Variance  power  spectra  of  wind  at  94  (upper)  and  106  (lower)  km 
during  5-15  August,  1976.  Solid  line  is  zonal  wind  and  dashed 
line  is  meridional  wind.  The  ordinate  is  percent  of  total 
variance,  which  is  given  on  the  upper  right  of  each  chart. 

Vertical  arrows  are  entered  at  the  point  where  50%  of  the  variance 
lies  beneath  the  curve  on  either  side  of  the  arrow. 
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Publications  and  Papers 


1.  The  following  paper  based  in  part  on  the  results  of  this  research, 
is  being  prepared  for  submission  to  Geophysical  Research  Letters: 
"The  QBO  in  thermospheric  wind  and  its  relation  to  the  QBO  in 
geomagnetism. " 

2.  An  extract  of  the  results  is  also  accepted  for  the  IAGA/IAMAP 
Joint  Assembly,  in  Seattle,  22  August  to  3 September  1977.  The 
title  is  " Long-period  waves  in  mesospheric  winds  at  Saskatoon 
(52°N) . " 


